\(\int \frac {1}{(d f+e f x)^2 (a+b (d+e x)^2+c (d+e x)^4)} \, dx\) [643]
Optimal result
Integrand size = 33, antiderivative size = 204 \[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=-\frac {1}{a e f^2 (d+e x)}-\frac {\sqrt {c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b-\sqrt {b^2-4 a c}} e f^2}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b+\sqrt {b^2-4 a c}} e f^2}
\]
[Out]
-1/a/e/f^2/(e*x+d)-1/2*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(1+b/(-4*a*c+b^2)^
(1/2))/a/e/f^2*2^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2)-1/2*arctan((e*x+d)*2^(1/2)*c^(1/2)/(b+(-4*a*c+b^2)^(1/2))^
(1/2))*c^(1/2)*(1-b/(-4*a*c+b^2)^(1/2))/a/e/f^2*2^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2)
Rubi [A] (verified)
Time = 0.22 (sec) , antiderivative size = 204, normalized size of antiderivative = 1.00, number of
steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.121, Rules used = {1156, 1137, 1180, 211}
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=-\frac {\sqrt {c} \left (\frac {b}{\sqrt {b^2-4 a c}}+1\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a e f^2 \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {\sqrt {b^2-4 a c}+b}}\right )}{\sqrt {2} a e f^2 \sqrt {\sqrt {b^2-4 a c}+b}}-\frac {1}{a e f^2 (d+e x)}
\]
[In]
Int[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x]
[Out]
-(1/(a*e*f^2*(d + e*x))) - (Sqrt[c]*(1 + b/Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt
[b^2 - 4*a*c]]])/(Sqrt[2]*a*Sqrt[b - Sqrt[b^2 - 4*a*c]]*e*f^2) - (Sqrt[c]*(1 - b/Sqrt[b^2 - 4*a*c])*ArcTan[(Sq
rt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a*Sqrt[b + Sqrt[b^2 - 4*a*c]]*e*f^2)
Rule 211
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]
Rule 1137
Int[((d_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_), x_Symbol] :> Simp[(d*x)^(m + 1)*((a + b*x^2 +
c*x^4)^(p + 1)/(a*d*(m + 1))), x] - Dist[1/(a*d^2*(m + 1)), Int[(d*x)^(m + 2)*(b*(m + 2*p + 3) + c*(m + 4*p +
5)*x^2)*(a + b*x^2 + c*x^4)^p, x], x] /; FreeQ[{a, b, c, d, p}, x] && NeQ[b^2 - 4*a*c, 0] && LtQ[m, -1] && In
tegerQ[2*p] && (IntegerQ[p] || IntegerQ[m])
Rule 1156
Int[(u_)^(m_.)*((a_.) + (b_.)*(v_)^2 + (c_.)*(v_)^4)^(p_.), x_Symbol] :> Dist[u^m/(Coefficient[v, x, 1]*v^m),
Subst[Int[x^m*(a + b*x^2 + c*x^(2*2))^p, x], x, v], x] /; FreeQ[{a, b, c, m, p}, x] && LinearPairQ[u, v, x]
Rule 1180
Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]
Rubi steps \begin{align*}
\text {integral}& = \frac {\text {Subst}\left (\int \frac {1}{x^2 \left (a+b x^2+c x^4\right )} \, dx,x,d+e x\right )}{e f^2} \\ & = -\frac {1}{a e f^2 (d+e x)}+\frac {\text {Subst}\left (\int \frac {-b-c x^2}{a+b x^2+c x^4} \, dx,x,d+e x\right )}{a e f^2} \\ & = -\frac {1}{a e f^2 (d+e x)}-\frac {\left (c \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{2 a e f^2}-\frac {\left (c \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right )\right ) \text {Subst}\left (\int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx,x,d+e x\right )}{2 a e f^2} \\ & = -\frac {1}{a e f^2 (d+e x)}-\frac {\sqrt {c} \left (1+\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b-\sqrt {b^2-4 a c}} e f^2}-\frac {\sqrt {c} \left (1-\frac {b}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a \sqrt {b+\sqrt {b^2-4 a c}} e f^2} \\
\end{align*}
Mathematica [A] (verified)
Time = 0.24 (sec) , antiderivative size = 209, normalized size of antiderivative = 1.02
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=-\frac {\frac {2}{d+e x}+\frac {\sqrt {2} \sqrt {c} \left (b+\sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {\sqrt {2} \sqrt {c} \left (-b+\sqrt {b^2-4 a c}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} (d+e x)}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {b+\sqrt {b^2-4 a c}}}}{2 a e f^2}
\]
[In]
Integrate[1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x]
[Out]
-1/2*(2/(d + e*x) + (Sqrt[2]*Sqrt[c]*(b + Sqrt[b^2 - 4*a*c])*ArcTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b - Sqrt[
b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b - Sqrt[b^2 - 4*a*c]]) + (Sqrt[2]*Sqrt[c]*(-b + Sqrt[b^2 - 4*a*c])*Ar
cTan[(Sqrt[2]*Sqrt[c]*(d + e*x))/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b + Sqrt[b^2 - 4*a*c]])
)/(a*e*f^2)
Maple [C] (verified)
Result contains higher order function than in optimal. Order 9 vs. order 3.
Time = 0.66 (sec) , antiderivative size = 172, normalized size of antiderivative = 0.84
| | |
| method | result | size |
| | |
| default |
\(\frac {\frac {\munderset {\textit {\_R} =\operatorname {RootOf}\left (c \,e^{4} \textit {\_Z}^{4}+4 c d \,e^{3} \textit {\_Z}^{3}+\left (6 c \,d^{2} e^{2}+b \,e^{2}\right ) \textit {\_Z}^{2}+\left (4 d^{3} e c +2 b d e \right ) \textit {\_Z} +d^{4} c +b \,d^{2}+a \right )}{\sum }\frac {\left (-\textit {\_R}^{2} c \,e^{2}-2 \textit {\_R} c d e -c \,d^{2}-b \right ) \ln \left (x -\textit {\_R} \right )}{2 e^{3} c \,\textit {\_R}^{3}+6 c d \,e^{2} \textit {\_R}^{2}+6 c \,d^{2} e \textit {\_R} +2 d^{3} c +b e \textit {\_R} +b d}}{2 a e}-\frac {1}{a e \left (e x +d \right )}}{f^{2}}\) |
\(172\) |
| risch |
\(-\frac {1}{a e \,f^{2} \left (e x +d \right )}+\frac {\left (\munderset {\textit {\_R} =\operatorname {RootOf}\left (\left (16 f^{8} e^{4} c^{2} a^{5}-8 b^{2} f^{8} e^{4} c \,a^{4}+b^{4} f^{8} e^{4} a^{3}\right ) \textit {\_Z}^{4}+\left (12 a^{2} b \,c^{2} e^{2} f^{4}-7 a \,b^{3} c \,e^{2} f^{4}+b^{5} e^{2} f^{4}\right ) \textit {\_Z}^{2}+c^{3}\right )}{\sum }\textit {\_R} \ln \left (\left (\left (40 a^{5} c^{2} e^{5} f^{8}-22 a^{4} b^{2} c \,e^{5} f^{8}+3 a^{3} b^{4} e^{5} f^{8}\right ) \textit {\_R}^{4}+\left (25 a^{2} b \,c^{2} e^{3} f^{4}-14 a \,b^{3} c \,e^{3} f^{4}+2 b^{5} e^{3} f^{4}\right ) \textit {\_R}^{2}+2 c^{3} e \right ) x +\left (40 a^{5} c^{2} d \,e^{4} f^{8}-22 a^{4} b^{2} c d \,e^{4} f^{8}+3 a^{3} b^{4} d \,e^{4} f^{8}\right ) \textit {\_R}^{4}+\left (4 a^{4} c^{2} e^{3} f^{6}-5 a^{3} b^{2} c \,e^{3} f^{6}+a^{2} b^{4} e^{3} f^{6}\right ) \textit {\_R}^{3}+\left (25 a^{2} b \,c^{2} d \,e^{2} f^{4}-14 a \,b^{3} c d \,e^{2} f^{4}+2 b^{5} d \,e^{2} f^{4}\right ) \textit {\_R}^{2}+2 c^{3} d \right )\right )}{2}\) |
\(376\) |
| | |
|
|
|
|
[In]
int(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x,method=_RETURNVERBOSE)
[Out]
1/f^2*(1/2/a/e*sum((-_R^2*c*e^2-2*_R*c*d*e-c*d^2-b)/(2*_R^3*c*e^3+6*_R^2*c*d*e^2+6*_R*c*d^2*e+2*c*d^3+_R*b*e+b
*d)*ln(x-_R),_R=RootOf(c*e^4*_Z^4+4*c*d*e^3*_Z^3+(6*c*d^2*e^2+b*e^2)*_Z^2+(4*c*d^3*e+2*b*d*e)*_Z+d^4*c+b*d^2+a
))-1/a/e/(e*x+d))
Fricas [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 1477 vs. \(2 (167) = 334\).
Time = 0.31 (sec) , antiderivative size = 1477, normalized size of antiderivative = 7.24
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=\text {Too large to display}
\]
[In]
integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="fricas")
[Out]
1/2*(sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/(
(a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2
*(b^2*c^2 - a*c^3)*d + sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)
/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sq
rt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4)))
- sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a
^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(
b^2*c^2 - a*c^3)*d - sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/(
(a^6*b^2 - 4*a^7*c)*e^4*f^8)) - (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(-((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt
((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) + b^3 - 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) -
sqrt(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*
b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(b^2
*c^2 - a*c^3)*d + sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^
6*b^2 - 4*a^7*c)*e^4*f^8)) + (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^
4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) + sqr
t(1/2)*(a*e^2*f^2*x + a*d*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^2
- 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))*log(-2*(b^2*c^2 - a*c^3)*e*x - 2*(b^2*c^2
- a*c^3)*d - sqrt(1/2)*((a^3*b^4 - 6*a^4*b^2*c + 8*a^5*c^2)*e^3*f^6*sqrt((b^4 - 2*a*b^2*c + a^2*c^2)/((a^6*b^
2 - 4*a^7*c)*e^4*f^8)) + (b^5 - 5*a*b^3*c + 4*a^2*b*c^2)*e*f^2)*sqrt(((a^3*b^2 - 4*a^4*c)*e^2*f^4*sqrt((b^4 -
2*a*b^2*c + a^2*c^2)/((a^6*b^2 - 4*a^7*c)*e^4*f^8)) - b^3 + 3*a*b*c)/((a^3*b^2 - 4*a^4*c)*e^2*f^4))) - 2)/(a*e
^2*f^2*x + a*d*e*f^2)
Sympy [A] (verification not implemented)
Time = 3.20 (sec) , antiderivative size = 258, normalized size of antiderivative = 1.26
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=\operatorname {RootSum} {\left (t^{4} \cdot \left (256 a^{5} c^{2} e^{4} f^{8} - 128 a^{4} b^{2} c e^{4} f^{8} + 16 a^{3} b^{4} e^{4} f^{8}\right ) + t^{2} \cdot \left (48 a^{2} b c^{2} e^{2} f^{4} - 28 a b^{3} c e^{2} f^{4} + 4 b^{5} e^{2} f^{4}\right ) + c^{3}, \left ( t \mapsto t \log {\left (x + \frac {- 64 t^{3} a^{5} c^{2} e^{3} f^{6} + 48 t^{3} a^{4} b^{2} c e^{3} f^{6} - 8 t^{3} a^{3} b^{4} e^{3} f^{6} - 10 t a^{2} b c^{2} e f^{2} + 10 t a b^{3} c e f^{2} - 2 t b^{5} e f^{2} + a c^{3} d - b^{2} c^{2} d}{a c^{3} e - b^{2} c^{2} e} \right )} \right )\right )} - \frac {1}{a d e f^{2} + a e^{2} f^{2} x}
\]
[In]
integrate(1/(e*f*x+d*f)**2/(a+b*(e*x+d)**2+c*(e*x+d)**4),x)
[Out]
RootSum(_t**4*(256*a**5*c**2*e**4*f**8 - 128*a**4*b**2*c*e**4*f**8 + 16*a**3*b**4*e**4*f**8) + _t**2*(48*a**2*
b*c**2*e**2*f**4 - 28*a*b**3*c*e**2*f**4 + 4*b**5*e**2*f**4) + c**3, Lambda(_t, _t*log(x + (-64*_t**3*a**5*c**
2*e**3*f**6 + 48*_t**3*a**4*b**2*c*e**3*f**6 - 8*_t**3*a**3*b**4*e**3*f**6 - 10*_t*a**2*b*c**2*e*f**2 + 10*_t*
a*b**3*c*e*f**2 - 2*_t*b**5*e*f**2 + a*c**3*d - b**2*c**2*d)/(a*c**3*e - b**2*c**2*e)))) - 1/(a*d*e*f**2 + a*e
**2*f**2*x)
Maxima [F]
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=\int { \frac {1}{{\left ({\left (e x + d\right )}^{4} c + {\left (e x + d\right )}^{2} b + a\right )} {\left (e f x + d f\right )}^{2}} \,d x }
\]
[In]
integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="maxima")
[Out]
-1/(a*e^2*f^2*x + a*d*e*f^2) - integrate((c*e^2*x^2 + 2*c*d*e*x + c*d^2 + b)/(c*e^4*x^4 + 4*c*d*e^3*x^3 + c*d^
4 + (6*c*d^2 + b)*e^2*x^2 + b*d^2 + 2*(2*c*d^3 + b*d)*e*x + a), x)/(a*f^2)
Giac [F(-2)]
Exception generated. \[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=\text {Exception raised: TypeError}
\]
[In]
integrate(1/(e*f*x+d*f)^2/(a+b*(e*x+d)^2+c*(e*x+d)^4),x, algorithm="giac")
[Out]
Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:Error index.cc index_gcd Error: Bad Argument ValueError index.cc index_gcd Error: Bad Argument ValueDone
Mupad [B] (verification not implemented)
Time = 9.83 (sec) , antiderivative size = 4339, normalized size of antiderivative = 21.27
\[
\int \frac {1}{(d f+e f x)^2 \left (a+b (d+e x)^2+c (d+e x)^4\right )} \, dx=\text {Too large to display}
\]
[In]
int(1/((d*f + e*f*x)^2*(a + b*(d + e*x)^2 + c*(d + e*x)^4)),x)
[Out]
- atan(((-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12
*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^
13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(
a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b*c^3*e^
12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^
11*f^6)*1i + (-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3
*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2
*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))
/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c
^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3
*d*e^11*f^6)*1i)/((-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1
/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^
2*c^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^
3*c^2*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(
1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^
5*b*c^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^(1
/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^
2*c^3*d*e^11*f^6) - (-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^
(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*
b^2*c^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*
b^3*c^2*d*e^13*f^10)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)
^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*
a^5*b*c^3*e^12*f^8)*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c - a*c*(-(4*a*c - b^2)^3)^
(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*
b^2*c^3*d*e^11*f^6) + 2*a^3*c^4*e^10*f^4))*(-(b^5 + b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c -
a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*2i - ata
n(((-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^
4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12*f^6)
- ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^13*f^
10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b
^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b*c^3*e^12*f^
8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^
4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^11*f^
6)*1i + (-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3*e^12
*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2*d*e^
13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(
a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c^3*e^
12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3*d*e^
11*f^6)*1i)/((-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c^3
*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c^2
*d*e^13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))
/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*b^3*c^2*e^12*f^8 - 16*a^5*b*c
^3*e^12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2))/
(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c^3
*d*e^11*f^6) - (-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2)
)/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*(x*(4*a^4*c^4*e^12*f^6 - 2*a^3*b^2*c
^3*e^12*f^6) - ((x*(8*a^5*b^3*c^2*e^14*f^10 - 32*a^6*b*c^3*e^14*f^10) - 32*a^6*b*c^3*d*e^13*f^10 + 8*a^5*b^3*c
^2*d*e^13*f^10)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2
))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) - 4*a^4*b^3*c^2*e^12*f^8 + 16*a^5*b
*c^3*e^12*f^8)*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(-(4*a*c - b^2)^3)^(1/2)
)/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2) + 4*a^4*c^4*d*e^11*f^6 - 2*a^3*b^2*c
^3*d*e^11*f^6) + 2*a^3*c^4*e^10*f^4))*(-(b^5 - b^2*(-(4*a*c - b^2)^3)^(1/2) + 12*a^2*b*c^2 - 7*a*b^3*c + a*c*(
-(4*a*c - b^2)^3)^(1/2))/(8*(a^3*b^4*e^2*f^4 + 16*a^5*c^2*e^2*f^4 - 8*a^4*b^2*c*e^2*f^4)))^(1/2)*2i - 1/(a*e*(
d*f^2 + e*f^2*x))